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ABSTRACT 


S - state excitation cross sections for proton-hydrogen 
atom collisions are calculated by a non-adiabatic method. 


' In .this letter we consider excitation of hydrogen atoms according to 
reactions of the type. 

p + H(is) - p + H(ns) (l) 

using a generalization of the nonadiabatic theory of Temkin (l) . Let r and 
R be the position vectors of the electron and incident proton relative to the 
(stationary) target proton. Above a few keV we may safely apply the impact 
parameter method which allows us to take R = £ + vt where v is the velocity 
of the incident proton and t is the time chosen such that at t - o the protons 
have a minimum separation, £. The electronic wave function Y (r,t ) satisfies 
the time-dependent Schroedinger equation (in atomic units) 

»-V* +i k ” -r {2) 


We have retained the interproton potential so that at large separations the 
perturbation tends rapidly to zero* 

Since ve are here interested only in the excitation of s- states j i 
natural to approximate Y by 

Where is the normalised hydrogenic s- state with principal quantum 

* AS 

number n and binding energy e n . Substituting ( 3 ) in (2) and integrating over 
the angular variables gives 

where 

= ° j r < ] 

_ 1 _i. 

r r, j 


t > R 


( 5 ) 
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Equation (4) must be solved subject to the boundary conditions 

L_\m ^ \S } V\ = L\n\ 5"^^^ = ° (6) 

o 00 

and for a hydrogen atom initially in the is state we must also have 

L im (c $ ~ Ri s ^ 

The problem defined above for closely resembles the zeroth order 

problem of Temkin's nonadiabatic theory. As in Temkin's application, equation 
(4) has a clear physical interpretation. Initially the electron sees only the 
charge of the target proton. However, when the electron cloud is penetrated 
by the incident proton, the electron sees a doubly charged nucleus. Thus a 
temporary helium atom is formed which decays as the moving proton emerges from 
the electron cloud. 

The excitation emplitudes can be obtained directly from equation 

(3) 

ot K it^ = e' £,vt Juan Ar (8) 

-Jo 


or, by making use of equation (4) they can be written in the integral form 


1 00 


oo 


R Vis irt i«- (9) 

RUM 


The boundary conditions at infinity (6) and (7), are awkward to handle 
numerically and it was found more convenient to make the transformations 


__ — 
5 = e 


X, t - US' (^ ) (10) 

thus placing the entire problem within a box 4 4 O ^ ^ I- 
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The boimdary conditions on x ma y now be simply stated 


X U;t') - ) 


(u) 


X (0, * ) = X If , \ } = 0 02) 

The finite (central) difference equation corresponding to (4) is 


A X = K 


(13) 


where A is a tri-diagonal complex matrix, y is a column vector representing 
the solution at some value of T and K depends on previous values of t. Using 
a non-itterative technique developed by one of us (E. S.) (2) equation (13) 
may be directly inverted and obtained from a knowledge of 

X (? y') • Since the initial condition at /v £= - -£• is given by 

(ll) we can develope a numerical solution over the entire region — 

Calculations were performed on the Laboratory's IBM 7094/7040 using a 
400 x 400 point mesh on the \ ^ plane. The unitarity requirement 


f /j - 


\ 1 % <*4 

Jo 


(14) 


was confirmed to a very high degree of accuracy except for values of t close 
to n/2 where it was occasionally violated by as much as 2 $. This was due to 
an accumulated instability in X which greatly inhibited the convergence of the 
excitation probabilities when computed using equation (8). 

However, the instability had an insignificant effect on the values of ^ 
computed from ( 9 ) which depends on the entire history of X rather than upon 
its instantaneous value . 
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The excitation cross sections, computed from 

CO 

Q 2j f d f <W> 

using the ^^0 obtained from (9) are given in the table, 

corresponds to the s-wave contribution to the ionization cross section and 
was computed from 

J 0 

where the complete sum was formed by extrapolating the computed results 
above n:7. 


1. A. Temkin, Phys* Rev* 126, 130 (1962) 

2* A. Temkin and E. C. Sullivan, NASA Technical Note D-1702 (unpublished) 
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Excitation 

Cross 

Sections in Units 

of Ara* 

/n 3 ) 

25 

50 

100 

200 

4oo 

800 

1.68 

1.17 

• 716 

.395 

.209 

.108 

1.22 

.827 

• 521 

.269 

.141 

.0725 

1.10 

.741 

.445 

.238 

.125 

.0639 

1.05 

.705 

.412 

.225 

.118 

.0604 

1.05 

.687 

.402 

.219 

.114 

.0587 

1.01 

.6 76 

.594 

.215 

.113 

.0576 

.269 

.222 

.135 

.0718 

0370 

.0188 


